The Absence of Vortex Lattice Melting in a Conventional Superconductor 
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The state of the vortex lattice extremely close to the superconducting to normal transition in an 
applied magnetic field is investigated in high purity niobium. We observe that thermal fluctuations 
of the order parameter broaden the superconducting to normal transition into a crossover but no 
sign of a first order vortex lattice melting transition is detected in measurements of the heat capacity 
or the small angle neutron scattering (SANS) intensity. Direct observation of the vortices via SANS 
always finds a well ordered vortex lattice. The fluctuation broadening is considered in terms of the 
Lowest Landau Level theory of critical fluctuations and scaling is found to occur over a large H C 2 (T) 
range. 



I. INTRODUCTION 



Thermal fluctuations are large in the high- T c super- 
conductors and their signature can be seen over a signif- 
icant region of the superconducting phase diagram. In 
high quality samples, the effects of disorder can be neg- 
ligible compared with thermal fluctuations and a clear 
first order transition, associated with melting of the vor- 
tex lattice 1 -! 2 - occurs well below the "mean-field" super- 
conducting to normal transition temperature, T C 2{H). 
However, in a low-T c sample, thermal fluctuations are 
significantly smaller due to the low temperature, the high 
vortex rigidity and strong inter-vortex forces. Using the 
Lindemann criterion (that the lattice melts when fluc- 
tuations reach a critical fraction of the lattice spacing) 
vortex lattice melting is predicted 3 - to occur only several 
millikelvin below T C 2{H). 

The small thermal energy in low- T c superconduc- 
tors makes pinning of the vortex lattice much more 
significant 4 , hence phase transitions related to this pin- 
ning become an alternative theoretical scenario^. For 
example, a disordering transition in the vortex lattice of 
niobium has been observedl directly by small angle neu- 
tron scattering (SANS) and is associated with the peak 
effect^ (increased pinning of the vortices, linked to a soft- 
ening of the vortex lattice below T '^(i?)). The position 
of this transition varies with niobium sample purity^ and 
is not observed, to within 20 mK of T C 2(H), in a high pu- 
rity sample 1 ^. Heat capacity measurements on NbaSn 11 
show signs of a vortex lattice phase transition in a sample 
with significant vortex lattice pinning (and showing the 
peak effect), hence it is not clear whether this transition 
is thermal melting or pinning-influenced disordering. To 
investigate thermal melting of the vortex lattice, a high- 
quality low-pinning sample is clearly important. Here we 
present heat capacity and SANS measurements on high 
purity niobium. We observe the effects of thermal fluctu- 
ations of the superconducting order parameter, but find 
no evidence for either a first order melting transition or 
for order-disorder transitions associated with vortex pin- 
ning. 



II. METHODS 

Heat capacity measurements were performed using the 
a.c. technique of Sullivan and Seidell. The sample 
is strongly thermally connected to a thermometer and 
heater, and weakly connected to a bath. An a.c. current 
of frequency \io is passed through the heater, so that 
heat is generated at a rate ^Pq [cos(u;t) + 1]. Conditions 
are such that: a) The heat capacities of the thermome- 
ter and heater are much less than that of the sample, 
b) The sample heater and thermometer come into equi- 
librium with a time constant much less than the inverse 
of the angular frequency, uj. c) The sample to bath re- 
laxation time is much longer than the inverse frequency. 
Under these conditions the temperature of the sample is 
given by a d.c. increase above the temperature of the 
bath, plus an a.c. component, 

T a .c. = P COs(Ljt)/(20jC), 

where C is the heat capacity of the sample (plus negli- 
gible contributions from heater and thermometer). The 
heat capacity can then be measured by accurate measure- 
ment of the size of the a.c. temperature oscillation. One 
important caveat of using this technique is that changes 
in the heat capacity over a range equal to the amplitude 
of the temperature oscillation are smoothed over. It is 
therefore important to set the size of the temperature 
oscillation to be small enough that sufficiently detailed 
changes in the heat capacity can be resolved. The size 
of the a.c. temperature oscillation at the sample was 
therefore set at less than 1 mK (typically 0.2 mK). 

A schematic of the cryostat insert is shown in Figure [1] 
A 14 g niobium sample (Nbl) was attached to a copper 
holder by nylon thread, and thermally connected to it by 
copper wire to give a thermal time constant of ^ 3 s to 
the holder. Heat was applied to the sample by passing 
a current at 4 Hz along a small amount of resistive wire 
wound non-inductively around the sample. Measurement 
of the oscillatory part of the sample temperature was 
performed using a ground-down GaAlAs diode; the a.c. 
component of the voltage across it was measured using 
a lock-in-amplifier (LIA) at 8 Hz. The sample holder 
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Figure 1: (color online) Schematic of the cryostat insert for 
the heat capacity measurements showing how the sample 
holder connects to the 1 K Pot via the temperature control 
plate. All of these were mounted inside a vacuum can, with 
no exchange gas. Heaters and thermometers were present on 
the sample, on the temperature control plate and on the 1 K 
Pot. 



consisted of a large block of copper to act as a low-pass 
thermal filter. The sample holder was strongly thermally 
connected to a temperature control plate, where the base 
temperature of the sample was controlled above that of 
the 1 K Pot during a measurement scan. The required 
temperature stability (~ 20 /iK) was achieved using a 
home-made analog PID (proportional, integral, deriva- 
tive) feedback system that heated the plate. The sensor 
was a high-sensitivity carbon resistance thermometer at- 
tached to the temperature control plate, with its resis- 
tance measured by an audio-frequency a.c. bridge with a 
LIA as null detector. Fine temperature scans were made 
by adding a ramped d.c. voltage to the LIA d.c. output 
to shift the control point of the PID controller. Abso- 
lute measurement of the temperature at the temperature 
control plate was achieved by a calibrated Ge thermome- 
ter mounted onto the plate (but away from the applied 
magnetic field). The sample temperature was corrected 
to take into account the d.c. component of the heating 
at the sample. 

During a measurement, the temperature was ramped 
with very small steps through the superconducting tran- 
sition at fixed field both on warming and on cooling. Nor- 
mal state measurements were taken 20 mT higher in field. 
A typical scan took 24 hours. A measurement of the am- 
plitude of the a.c. voltage from the GaAlAs thermometer 
was typically taken every 300 s, with the lock-in-amplifier 
averaging the signal with a time constant of typically 100 
s. This averaging introduced a lag between warming and 
cooling measurements which has been corrected for in the 
measurements presented in the paper. 

The sample was aligned in the centre of a magnetic 
field with uniformity 3 parts in 10 4 over the volume of the 
sample. Two small superconducting coils were arranged 
either side of the sample to allow an optional a.c. field of 
amplitude 0.7 G to be applied perpendicular to the main 
field. 

Small angle neutron scattering (SANS) measurements 
were performed on the instrument D22 at the Institut 
Laue Langevin, Grenoble but with the instrument tem- 



perature control enhanced to achieve the required ultra- 
fine temperature control. While the external electronics 
were the same as the heat capacity measurements, the 
nature of neutron cryostats required a different set-up in 
the cryostat insert. There was a low pressure of helium 
exchange gas in the sample chamber and its tempera- 
ture was set by the relatively coarse temperature control 
of the instrument. The sample was mounted in a cop- 
per can inside the sample chamber; the can temperature 
was raised above that of the instrument control using a 
heater and carbon resistance thermometer attached in- 
side the can, and PID feedback was used, as in the heat 
capacity measurements to perform temperature sweeps. 
The magnet that provided the main magnetic field was 
designed to give uniformity 1.5 x 10 -4 over a 24 mm 
diameter spherical volume. As with the heat capacity 
measurements, additional coils allowed a small a.c. field 
of amplitude 0.1 G to be applied to the sample. 

Our experiments required samples in which vortex lat- 
tice pinning was minimal. The 14 g, barrel shaped, nio- 
bium sample (residual resistance ratio, RRR = 1000) 
was made high-purity and homogeneous by annealing in 
ultra-high vacuum (< 3 x 10~ 9 torr) for 48 hr at 2300 K. 
It was then oxygen treated^ to reduce surface pinning. 
This sample was used in both the heat capacity measure- 
ments and SANS measurements, and is known as Nbl. 
Two other samples (used in the SANS measurements), 
Nb2 and Nb3, were spherical, 11 g (RRR = 450) and 
ellipsoidal, 8 g (RRR = 1000) respectively, and prepared 
in the same way. 



III. HEAT CAPACITY MEASUREMENTS 

Measurements of the heat capacity of niobium were the 
first to show that critical fluctuations broaden T C 2(H) 
from a second order phase transition into a crossover—. 
While no first order vortex lattice melting transition was 
observed in these measurements, they were performed at 
constant magnetic induction B, whereas constant mag- 
netic field H is required to observe a first order transi- 
tion clearly. Figure [5] presents our measurements of the 
heat capacity of high purity niobium over the full H C 2(T) 
range. The surface treatment given helps to make the 
sample almost perfectly magnetically reversible, to al- 
low the measurements to be performed at constant H. 
However, a small (< 0.5 mK) hysteresis can be seen at 
intermediate fields, resulting from some residual surface 
pinning of vortices and corresponding to flux trapping 
of less than 0.1 mT. Magnetisation measurements on a 
smaller sample, similarly treated, showed no sign of the 
peak effect. 

As shown in Figure [3J we observe a transition broad- 
ened by critical fluctuations over the range ~ 5 to 30 mK 
(depending on field), giving a peak below the mean field 
transition and a long tail above it. At no temperature do 
we see sign of a first order phase transition, of the type 
observed in the high- T c superconductors, which would 
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Figure 2: (color online) The ratio of the superconducting and normal state heat capacities, C a /C„, measured as a function 
of temperature in the niobium sample Nbl. The magnetic field was applied parallel to the [111] crystal axis. Normal state 
measurements, C n , were taken 20 mT higher in field. The temperature has been plotted as T-T C 2,mft, where T C 2,mft(H) is 
the value the transition temperature would have in the absence of fluctuations and is taken as the temperature at which the 
fluctuation heat capacity is a fraction, 0.2, of the mean-field jump in C~. When a small hysteresis is seen, measurements taken 
on warming show the higher transition temperature. 



appear as a spike in the heat capacity. In YBa2Cu307_,5 
the entropy of vortex lattice melting was found to be 0.4 
ks per vortex per superconducting layer 2 . In comparing 
the latent heat in a layered system with that in isotropic 
niobium, we take the layer spacing to be equivalent to 
the coherence length, £. Using parameters suitable for 
our niobium sample we estimate (Appendix) that melt- 
ing would cause a spike in the heat capacity of height 
AC meit - 1 mJ mol" 1 K" 1 in a field of 300 mT (for 
a transition width AT me i t ~ 2 mK); this is ~ 10 % of 
the mean field jump, AC, at the superconducting to nor- 
mal transition, which would be readily resolvable in our 
measurements but is clearly not seen. 



IV. SANS MEASUREMENTS 

We can also search for a melting transition by examin- 
ing the state of the vortex lattice using SANS. This tech- 
nique was used in the observation of melting/ decoupling 
of 'pancake' vortices in the high-T c material BSCCO 15 
but the low scattering intensity allowed no investiga- 
tion of the nature of the disordered vortex state. The 
straightforward geometry of SANS allows the size of the 
diffraction spot on the detector to be related to corre- 
lation lengths of the vortex lattice 16 . This idea is illus- 



a "rocking curve" of diffraction intensity as a function of 
angle when </>, the angle between the applied magnetic 
field/sample [111] axis and the incident neutron beam, 
is rocked through the Bragg condition. On the peak of 
the rocking curve, uj rac i and cutang are the widths of the 
diffraction spot on the detector, parallel and perpendicu- 
lar respectively to the relevant g-vector. The collimation 
and wavelength spread contribute to these widths, but 
there are additional contributions due to the crystallinity 
of the vortex lattice in the sample. Mosaic spread (i.e. 
meandering of the average vortex line direction) is mainly 
seen in Ui n . Any spread in lattice spacing is more appar- 
ent in ujrad, while the spread in vortex lattice orientation 
(corresponding to rotations of the vortex lattice about 
the field direction) contributes to uJtang- 

Measurements of these widths, on the same sample 
used in the heat capacity measurements, Nbl, are pre- 
sented in Figure^b). The perfection of the vortex lattice 
is such that the width of the rocking curve, uJi„, is lim- 
ited by the instrument resolution while no changes in the 
radial spot width Lu ra d are observed. There is a small 
increase in cutang, the tangential width of the spot, but 
the ~ 0.01° change is extremely small compared with 
the spread that would be observed for an orientation- 
ally disordered vortex lattice, which would give a ring of 
diffracted intensity. 

trated in Figure [3£a) . The width uii„ is obtained from A signature of melting would also appear in the inten- 
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Figure 3: (color online) (a) Schematic to illustrate the geometry of a SANS experiment, where ki and k/ represent the incident 
and scattered neutron beams respectively, (b) The FWHM of Gaussian fits to the size of diffraction spots at the detector in the 
temperature regime of critical fluctuations at 200 mT (applied parallel to the [111] axis of Nbl). Measurements were taken on 
cooling and are shown for all temperatures at which there was sufficient diffracted intensity to fit the spot shape. Backgrounds 
taken at 5.610 K were subtracted before fitting. For measurements of uj ra d and uitang, <t> was fixed at the angle that brought the 
top two first order spots from the hexagonal lattice, left (#) and right (*), simultaneously onto the Bragg diffraction condition 
and maximised their intensity. Red lines mark the FWHM expected for a perfect lattice, where the width is due solely to 
instrument resolution. 



sity, Ihk, of a diffraction spot, integrated over the rocking 
curve. This depends on the vortex lattice form factor, 
Fhk, which expresses the "field contrast" in the mixed 
state, 

where /i is the incident neutron flux, V is the sam- 
ple volume, 7 is the neutron magnetic moment in nu- 
clear magnetons, A n is the neutron wavelength and $o is 
the flux quantum. Mean-field Ginzburg-Landau theory^ 
gives Fhk proportional to the magnetisation, M, 

F hk = (-lf e -™/^VoAf, 

where v = h 2 + hk + k 2 . In the absence of fluctuations 
the magnetisation, and hence the form factor, should go 
linearly to zero as T — > T c2 (H). This linear variation of 
form factor with temperature has been verified to within 
20 mK of T C 2(H) on a high purity niobium sample 10 , 
however that experiment did not have sufficient temper- 
ature resolution to probe the region where vortex lattice 
melting would be expected to occur. Measurement of the 
integrated intensity this close to T C 2(H) is only possible 
due to the relatively large field contrast in the mixed 
state of niobium. The equivalent experiment on a super- 
conductor with a longer penetration depth would require 
unfeasibly long count times^. 

As shown in Figure 2(b), the rocking curve width, u>i n , 
for the vortex lattice in our sample is essentially indepen- 
dent of temperature in the critical region. Hence, apart 



from a proportionality constant, the temperature varia- 
tion of \Fhk\ may be obtained by tilting the sample to 
the peak of the rocking curve and measuring diffraction 
intensity versus temperature with <f) fixed. In the critical 
fluctuation region, we observe (Figure@]) that the square- 
root of the diffracted intensity, Vi, falls below the linear 
mean-field theory prediction, but even in this region the 
lattice remains well ordered (inset Figure[4]). A first order 
melting transition would appear as a discontinuity in the 
intensity, and this is not observed. It might be argued 
that temperature, field or sample inhomogeneity could 
provide the smearing of the SANS results. However, if 
we suppose that the measurements represent a smearing 
of mean field linear variation of \fl with temperature 
near T C 2 the intensity would lie above the mean field ex- 
trapolation from lower temperatures, and not below it as 
we observe. 

V. CRITICAL SCALING 

To examine the origin of these deviations from mean 
field theory we look for a universal temperature scale 
which would allow measurements at different fields to 
collapse onto the same curve. 3D XY scaling is not 
suitable, as it requires a spatially uniform order param- 
eter, which is not valid in niobium as the vortex cores 
are large. However, the Lowest-Landau-Level Ginzburg- 
Landau (LLL-GL) theory of critical fluctuations is appro- 
priate. The relevant temperature scale ax, as calculated 
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Figure 4: (color online) The square root of the diffracted intensity minus background, \/7 ( oc \Fhk\)j as a function of temperature 
in a field of 200 mT, measured during the same sca.ii that co r ad and cotang (Figure 2(b)) were recorded. Negative values of vT 
(calculated as sign(Z) yj/j) result from the background subtraction. The mean field theory linear behaviour is obtained from 
a fit to the temperature dependence outside the critical region. Insets show the diffraction pattern observed on the detector 
(with direct beam masked and background subtracted) for two temperature points. 



by Thoulessi^, is given by or = (T — mft)/S, where 
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Critical scaling of the heat capacity measurements, Fig- 
ure [5]^ a), is successful over a large region of the phase 
diagram, except at the highest and lowest fields. This 
confirms the fluctuation origin of the signal and rules 
out inhomogeneity broadening. In the field region where 
our measurements scale, we find that the long tail in heat 
capacity at high temperatures agrees well with the Thou- 
less prediction within LLL-GL, but the enhanced peak is 
narrower than expected from theory. This may demon- 
strate the need for higher-order terms in the fluctuations 
to be included in the theory, mixing of Landau levels at 
the transition^, or a consequence of coupling between 
the vortex and crystal lattice (which also breaks the ro- 
tational degeneracy of the vortex lattice and gives it a 
preferred orientation) . 

As the neutron scattering intensity originates from the 
magnetisation, it would be of interest to see how its be- 
haviour compares to that of the heat capacity. The heat 
capacity and SANS measurements were taken using dif- 
ferent magnets, and different thermometers, neither of 
which had calibration accuracy - in absolute tempera- 
ture - small compared with the widths of the transitions. 
However, by comparison of the data with theory, we have 



been able to inter-compare data taken at different fields 
and temperatures - and data from the two different tech- 
niques - as a function of temperature difference from 
the mean field transition temperature, T c2t MFT, even 
though there is no sharp transition once fluctuations are 
taken into account. 

To calculate the magnetisation of the sample from the 
heat capacity, we use the thermodynamic relation 
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Integrating the Thouless theoretical prediction for the 
heat capacity^ down from high temperatures, we obtain 
Figure (Hla), which indicates that at temperatures be- 
low the critical region the magnetisation varies linearly 
with temperature, which is the mean field GL result in 
the mixed state. However, it will be noticed that the 
magnetisation extrapolates linearly to zero at a reduced 
temperature of ay = -1.0. This indicates, as one would 
expect, that fluctuations not only broaden the transition, 
but also suppress it below the mean field value. In Figure 
IH^b) , experimental heat capacity data is integrated in the 
same way and a similar graph is obtained. Finally, we 
consider the SANS data, where below the critical region 
the SANS intensity obeys the mean-field result \fl oc M. 
The SANS intensity measurements can be placed on the 
same temperature scale by taking the extrapolation of 



6 




Figure 5: (color online) Scaling of the heat capacity and SANS intensity, (a) The heat capacity measurements of Figure 1, 
normalised to the mean-field behaviour as (C s - C n ) / (Cs.mft - C„), plotted on the Thouless reduced temperature scale, &t- 
C ' s,mft was extrapolated through the transition from a linear fit of C a / C n versus temperature below the fluctuation region^. 
The measurement at 44 mT is shifted vertically by 0.2 units, to avoid overlap with the Thouless theoretical prediction—^, which 
is plotted as the black line. All successive plots (at the same sequence of fields as in Figure 1) are shifted by O.f units, (b) 
SANS intensity, normalised as / / I mft, where Imft was found from a linear fit of vT versus temperature below the critical 
region (as shown for example in Figure 2) plotted versus the Thouless reduced temperature &t- Measurements are shown for 
three samples, Nbl: barrel shaped, Nb2: spherical and Nb3: ellipsoidal. All samples were prepared under similar condition 
with similar RRR values. As no significant difference was observed between the left and right hand diffraction spots, or on 
warming and cooling, these measurements have been plotted with the same marker style for clarity. The error bars become 
large at high temperatures because the data points represent the ratio of two quantities, both of which are tending to zero. 



the mean field behaviour, where \fl — > at &t = -1.0, 
as shown in Figure |^c). 

We can now compare the temperature-dependence of 
the magnetisation obtained from these two different tech- 
niques. The magnetisation derived from the heat capac- 
ity measurements demonstrates superconducting fluctu- 
ations result in a long tail above T C 2,mft, giving in this 
region a diamagnetic magnetisation larger than the mean 
field extrapolation. However, the SANS intensity de- 
creases faster than mean field behaviour and is already 
immeasurably small below T C 2,mft- The different tem- 
perature dependencies is because the heat capacity mea- 
surements are sensitive to all fluctuations that give rise 
to an average magnetisation, while the SANS intensity 
is sensitive only to a spatially coherent contrast in the 
magnetisation. Fluctuations of the superconducting or- 
der parameter are expected to result in local field fluc- 
tuations in a 1ow-k superconductor such as niobium 18 . 
Hence these local field fluctuations reduce the spatially 
coherent magnetic field contrast, causing the decrease in 
SANS intensity below the mean-field value. This is con- 
firmed in Figure Efb) , where the SANS intensity in the 
critical region is shown to obey scaling once plotted on 
the &t temperature scale. 



VI. OSCILLATING FIELD 



We also performed the heat capacity and SANS mea- 
surements with a small oscillating magnetic field applied 
to the sample in addition to the main field. This type 
of field can help the vortices find their equilibrium po- 
sition in the case of weak pinning. Figure UJa) shows 
the raw lock-in-amplifier output from the heat capacity 
measurements (the signal is oc 1/C and the output has 
not been corrected for time lag due to LIA time con- 
stant) taken under identical conditions except for the 
presence/absence of a small oscillatory magnetic field. 
The size of the magnetic field was too small to give a 
noticeable change in transition temperature in this scan 
(and because it was applied perpendicular, the value of 
the field was essentially unaltered; even parallel, A£> C 2 
~ 0.7 G corresponds to AT C 2 ~ 0.2 mK). The measure- 
ments with and without the oscillatory field are identical. 
Figure [7(b) shows a similar experiment performed on the 
SANS measurements, except in this case the oscillatory 
field was applied between each measurement, and con- 
sisted of a decaying exponential oscillatory field, with 
initial amplitude 0.1 G. No difference in diffracted inten- 
sity between these two measurements can be observed 
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Figure 6: (color online) The magnetisation extracted from 
the heat capacity and the SANS intensity, (a) C s - C n from 
Thouless theory integrated with respect to temperature plot- 
ted versus the Thouless reduced temperature, ar. (b) C s - 
C n from experiment, integrated with respect to temperature 
and plotted versus &t- Measurements taken from Nbl in 198 
mT. (c) The SANS diffraction intensity, presented as %/7 (oc 
M) versus &t- Measurements taken from Nbl in 200 mT. 
Red lines show the extrapolation of linear mean-field-theory 
behaviour. 



above the noise. That we observed identical behaviour, 
both with and without this additional oscillating field, 
illustrates the very low pinning nature of our sample and 
confirms that all our measurements were performed on a 
vortex lattice in equilibrium. 



VII. CONCLUSIONS 

In summary, measurements of the heat capacity and 
SANS intensity in high purity niobium show no sign of 
a first order vortex lattice melting transition. Direct ob- 
servation of the vortices by SANS always finds a well or- 
dered lattice. Fluctuations are observed to broaden the 
second order superconducting to normal transition into a 
crossover in both the heat capacity and SANS intensity 
measurements. Lowest Landau Level (LLL) scaling is 
successful over a large range of fields, characterising the 
nature of the fluctuations. The fluctuation enhancement 
of the heat capacity is found to be more sharply peaked 
than the predictions of Thouless 18 . 

We now have to consider why the vortex lattice does 
not melt in niobium, a task made difficult due to the 
paucity of theoretical work. Our measurements show 
that the temperature range over which thermal fluctu- 
ations affect the superconducting order parameter is in 
agreement with LLL theory. Calculations within LLL 
theory predict vortex lattice melting to occur at ay = -9.5 
but in a strongly Type-II superconductor—. This is well 
below our fluctuation region. How this prediction is mod- 
ified for Iow-k niobium is a current theoretical challenge. 
One difference is that in the high-T c superconductors the 
melting temperature is well below the mean- field T C 2 (H) 
so the order parameter is not changing much in magni- 
tude with temperature, however, in niobium the melting 
temperature is expected to be very close to T C 2(H) and 
therefore the loss of phase coherence of the order param- 
eter would occur in a region where it is rapidly changing. 

The absence of melting may be due to interactions 
of the vortex lattice with the underlying crystal lattice, 
which is known to have a significant influence on vor- 
tex behaviour at lower temperatures. For example, frus- 
tration between the ideal hexagonal vortex lattice and 
the cubic crystal lattice results in square, scalene and 
isosceles vortex lattice structures when the field is applied 
along a [001] crystal axis of niobium. These structures 
have been confirmed to persist up to H C 2(T) to within 
100 mK of T C 2(-ff)2i. These strong interactions with the 
crystal could maintain vortex lattice rigidity, preventing 
melting. It would be interesting to expand our measure- 
ments to applied field directions other than parallel to 
the [111] crystal axis, to examine the importance of this 
interaction in our fluctuation temperature regime. 

A more unconventional explanation for the absence 
of vortex lattice melting in niobium is also a possibil- 
ity. While there is no doubt that a first order transition 
occurs in the high-T c superconductors, the microscopic 
details of the transition are still poorly understood, es- 




Figure 7: (color online). Investigation of a small oscillatory field applied in addition to the main field, (a) Raw LIA output 
(oc 1/C) for Nbl in a field of 308 mT. In one measurement an oscillating field of amplitude 0.7 G, frequency 100 Hz was 
continually applied to the sample during the measurement, (b) Diffracted intensity from Nb3 taken with and without a small 
decaying oscillatory magnetic field applied to the sample before each measurement. The main field applied was 300 mT. 



pecially as direct observation of the vortices above the 
transition has not been achieved. Our results are com- 
patible with the controversial proposal^ that the vortex 
lattice is actually a liquid, but with a correlation length 
so large that it is experimentally indistinguishable from a 
solid, and the melting transition observed in the high- T c 
superconductors is related to decoupling of pancake lay- 
ers and would not be expected to occur in our isotropic 
system. 



Acknowledgments 

This work was supported by the UK EPSRC. Neutron 
experiments were performed on instrument D22 at the 
Institut Laue Langevin, Grenoble, France. We are grate- 
ful for discussions with M. A. Moore, B. Rosenstein and 
B. Ya. Shapiro. 



Appendix 

Vortex lattice melting in YBa2Cu307_^ was accompa- 
nied by an entropy change of A5* = 0.4 k# per vortex 
per layer—. We wish to make an estimate of the entropy 
change that would accompany vortex lattice melting in 
niobium. In comparing the latent heat in a layered sys- 
tem (where the vortex layer spacing is taken to be the 
distance between the Cu02 planes) with that in isotropic 
niobium, we take the layer spacing to be equivalent to 
the coherence length, so in niobium we expect A5* = 0.4 
ks/vortex/£, with £ = 29 nm. (The work of Dodgson et 
ali^ suggests that the relevant length scale is the vortex 
spacing and not £; however their theory applies in the 
London limit, and in the case considered below there is 



only a factor ~ 2 difference between these two lengths). 
The entropy density, As, is related to the entropy per 
vortex per layer, AS, by AS = d 2 £As, where d is the 
inter- vortex spacing and £ is the coherence length/layer 
spacing. In a field of 300 mT (0.7H c2 (T = 0), giving 
T C 2= 3.5 K) the inter-vortex spacing is d = 7.7xl0 -8 
m. This gives an entropy density of As = 32 mJ m 
K~ 3 = 3.48 x 10~ 4 mJ moh 1 K -1 . 

We can make an estimate for the height of the heat ca- 
pacity spike by assuming the spike is triangular in shape, 
with a transition width of AT me ; f . The width of this 
hypothetical melting transition will be related to sample 
homogeneity. In zero field (where fluctuations are only 
expected to be present on the scale of nK) our T C 2(0) 
width of 0.4 mK (measured from the heat capacity) is 
testament to the high quality of the sample. In a mag- 
netic field we need to consider how much any variations 
in sample purity will broaden T C 2, (additional to any 
broadening due to fluctuations of the order parameter). 
The upper critical field, B C 2, depends on the Ginzburg- 
Landau parameter, K\, as 5 c2 (T) = \/2Kii? c (T), and k\ 
depends on the electronic mean free path as given 2 ^ by 
k = kq + 7500poa/7i where kq is for the case of infinite 
mean free path, po is the residual resistivity and 7 is the 
Sommerfeld constant. Any variations in sample purity 
will give a spread in residual resistivity across the sam- 
ple, resulting in a spread in upper critical field, AB C 2 
~ (£,Al/l 2 )B C 2, where I is the electron mean free path. 
This estimation emphasises the importance of using a 
high-purity sample. 

We can relate 2 ^ electron mean free path to the resid- 
ual resistivity ratio, R, using I = 5 X 10~ 9 i?, so for our 
sample with R = 1000, I = 200£. Therefore, using AT c2 
~ {^Al / l 2 )B r2 {dB c2 / AT)^ 1 , we have, for an assumed 
inhomogeneity in Al of 5 %, a transition width due to 
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inhomogeneity broadening of AT C 2 ~ 2 mK in a field 
of 300 mT. If the width of the melting transition is the 
same: AT me i t ~ 2 mK, we can therefore estimate the 
height of the spike in the heat capacity to be AC me ; t <~ 



1.2 mJ mor 1 K" 1 in a field of 300 mT. The size of the 
mean field heat capacity jump in this field is ACmft 
= 10 mJ mol -1 K _1 so the spike would be 12 % of the 
jump. 
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